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We present the estimates of the Cosmic Microwave Background (CMB) polarisation power 
spectrum from galaxy clusters and filaments using hydrodynamical simulations of large scale 
structure formation. We focus in the present study on the CMB Quadrupole induced polari- 
sation, which is the dominate secondary polarisation effect, and give the E and B mode power 
spectra between I ~ 560 to 20000. 



1 Introduction 



After recombination, the CMB photons scatter off free electrons in ionised matter, such as the 
hot gas trapped in the potential wells of galaxy clusters. The presence of the CMB temperature 
quadrupole induces a linear polarisation in the scattered radiation. 

Sunyaev and Zel'dovich'^ were the first to estimate the level of polarisation in galaxy clus- 
ters. These authors pointed out that in addition to the primary CMB quadrupole there are two 
other sources of temperature quadrupole seen by a cluster: a quadrupole due to the transverse 
peculiar velocity of the cluster and double scattering. Studying the induced CMB polarisation 
due to clusters can open a whole new window for cosmology. Sunyaev and Zel'dovich proposed 
to use the polarisation to estimate cluster's transverse velocity. Measuring the polarisation to- 
wards distant clusters should provide us with an opportunity to observe the evolution of the 
CMB quadrupole^. Moreover, CMB quadrupole seen by the clusters contains statistical infor- 
mation on the last scattering surface of at the cluster p ositi on. Therefore, measuring the cluster 
polarisation should help us to beat the cosmic variance'^'. 

In this work, we study the polarisation induced by galaxy clusters. We focus in the statistical 
properties of such a polarised signal by computing its angular power spectrum at small angular 



scales. We explore the signal associated with the primary CMB temperature quadrupole since 
it dominates at the scales we are interested in^l. To do so, we follow the method of Liu et al.^l, 
in which the authors computed the polarised signal at reionisation using N-body simulations in 
combination with an analytical description to model the gas distribution. Based on the same 
formulae and similar approach, we estimate the polarisation from the hot gas in galaxy clusters 
and the warm gas in filaments using hydro dynamical simulations of large scale structure (LSS), 
which directly account for the gas dynamics. 

In Section 2 we describe the method of this work. We briefly review the formulae of the 
polarisation power spectra and the hydro dynamical simulations used in this paper. In Section 
3, we present the results of our work and give a discussion and conclusion in Section 4. 

2 Background 

The polarisation signal is usually described by the Stokes parameters which can be combined to 
obtain a divergence free component, the so-called E mode, and a curl component, the so-called 
B mode. 

The power spectra of polarisation of E and B modes are obtained by integrating the product 
of the visibility function g and the polarisation source. In the case of primary polarisation, 
the source is the quadrupole of the primary CMB temperature anisotropies. This is a first 
order signal. When dealing with second order effects, one has to take into account possible 
modulations of the visibility function by the electron density perturbations. This was done in 
Liu et al.^. More specifically, the power spectra were found to be: 
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where 5 e (k, r) is the electron density perturbation as a function of scale k and conformal time 
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! (k) is the projected scalar mode of primary CMB 

quadrupole (A^ - "' 1 ) in the frame of 5 e . The visibility function g(r) has a simple physical 
meaning, being the probability that a photon had its last scattering at epoch r and reaches 
the observer at the present time, tq. It is defined as g(r) = — ^e K ( T °)~ K ( T \ where k(t) is the 

electron-scattering optical depth at time r. The term T^ B y(kr) is given by the combination 

of the spherical Bessel functions (for details see Table 1 in Liu et al.^j). 

In order to calculate the polarisation anisotropy spectra using the formulae mentioned above, 
we compute the visibility function and the electron density distribution from clusters and fil- 
aments using non-radiative hydrodynamical simulations of a ACDM cosmology, with matter 
density parameter, J7 m = 0.3, cosmological constant density, = 0.7, hubble parameter, 
h = 0.71, and baryon density, Qf, = 0.044. The initial density field was constructed using 160 3 
particles of both baryonic and dark matter, perturbed from a regular grid of fixed comoving 
size L = 100/i -1 Mpc. The CDM matter power spectrum used in the simulation had a shape 
parameter, T = 0.21, and normalization eg = 0.9. The run was started at r edshift z = 50 and 
evolved to z = using the publicly available Hydra (AP 3 M/SPH) codeU^. The gravitational 
softening was fixed at 25/i~ 1 kpc in physical units below z = 1, and above this redshift held 
constant to 50/i _1 kpc in comoving co-ordinates. 

In the simulation, the gas is assumed ionised if the temperature is higher than 10 5 K. We 
further define two gas phases based on the baryon collapsed density: when the gas overdensity, 
5, is greater than the density contrast at collapse, A c = 178S7 m (,z)~ a55 , the gas is in clusters; the 
low density phase defined by 5 < 5 < A c encompasses a warm gas which we call filaments. The 
visibility function is obtained from the time evolution of the ionisation fraction in our run for 



T > 10 5 K. For the computation of C^,B)h we evaluate the electron density fluctuations on a 
regular grid with 600 3 cells inside the comoving box, which corresponds to a fixed cell separation 
of 0.17/i _1 Mpc in comoving coordinates. 

Finally, the only remained unknown term is the quadrupole of the temperature anisotropies. 
Given the cosmological para met ers and the ionisation fraction from our simulation, the pub- 
licly available code CMBFast — can output the time evolution of the temperature quadrupole 
component with different k-basis. Note that the perturbations of tensor mode are such that 
the normalizations to COBE and to as = 0.9, used in our simulation, are consistent. We as- 
sume that the tensor spectral index satisfies the single field inflationary consistency condition 
nt = —r/8 with r being the ratio of tensor to scalar quadrupole r = Cf jCf . Here r = 0.05 and 
nt = —0.00625 are applied when running CMBFast. 

3 Power spectra of polarisation of clusters and filaments 

We now investigate the contribution to the polarisation power spectra from the two different gas 
phases defined previously. We show the results in FigQ. Note that the range of multipoles i is 
limited by the resolution (for high is) and box size (for low is) of the simulation. The highest 
multipole i is set by the corresponding angular size of 0.17/i _1 Mpcat z = 0.72 in our simulation. 

We get the same power for the secondary E and B mode polarisation with a maximum 
difference is less than 10~ 6 . The reason for the equality of the total E- and B-mode power 
spectra is essentially that the first-order quadrupole whose scattering produces the polarisation is 
dominated by large scales, and so has a random orientation relative to the small-scale fluctuations 
in the electron density. Scattering of the quadrupole by the small-scale fluctuations therefore 
on average excites E- and B- modes equally 1 ^'. 

We now compare the polarisation signal from two phases: the cluster and so-called filaments. 
The power spectrum from all the ionized particles (with T > 10 5 K) is also shown in the figure 
(|T|). The two curves for all ionized particles and filaments have similar shapes peaking at 
I ~ 10 4 . The power spectrum of polarisation in clusters also peaks at the same scales but shows 
much less power on the large scales. The power spectrum of the filaments is about one order 
of magnitude larger than that of the clusters at I ~ 1000. This is due to the fact that the 
dominant contributions comes from the filaments on large scales. At smaller angular scales, 
the power spectrum of filaments is still larger than that of clusters. Recall the equation^, the 
power spectra is the integral of electron density weighted by the visibility function. From our 
simulation, the visibility function peaks at z >3 and its full width half maximum is bounded 
by z ~ 1.2 and z ~ 6. In this range, the power of electron density for filaments is quite larger 
than for clusters. In the other words, most of the scattering producing the linear polarisation 
occurs for z >1.2 when the relative contribution of the clusters is still small (for details see Liu 

et al.0). 

4 Discussion and Conclusion 

Using hydrodynamical LSS simulations, we compute the E and B mode polarisation power 
spectra of galaxy clusters and filaments due to the primary CMB quadrupole. As expected, 
these power spectra dominate at small angular scales. The question is then: will the secondary 
polarisation signal, associated with clusters and filaments, be a problem for the measurement of 
the primary E and B mode power spectra? To compare the polarisation spectra from the differ- 
ent physical processes at different scales, we show, in addition to our computations, the primary 
E and B modes by scalar and tensor mode perturbations respectively. We also show in the 
same figure the expected B mode signal from the power conversion of E mode by gravitational 
lensing. In our results, the E mode signal from clusters and filaments dominates the primary 
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Figure 1: Polarisation power spectra for all ionized particles, filaments and cluster are shown by dotted, dot- 
dashed and short dashed curves, respectively. We also plot the primary E mode induced by scalar perturbation 
(solid line), primary B mode (triple dot-dashed line) from tensor mode and B mode from lensed E mode (long 

dashed line). 

signal at / > 6000. Below this scale, the secondary signal is always very small and therefore is 
not expected to significantly contaminate the primary plarisation signal. 

For the B mode polarisation, and although some uncertainties in amplitude of the tensor 
mode perturbations, our secondary B mode signal is definitively not a problem for finding the 
gravitational wave induced polarisation at scales I <100. However, the gravitational lensing- 

induced signal which peaks around I ~ 1000 is a major issue^l The power spectrum of the 
B mode polarisation from clusters and filaments is two orders of magnitude smaller than the 
lensing and it dominates at smaller angular scales. It might be observed at these scales if the 
lensing-induced B modes can be significantly cleaned, see for example Seljak U. and Hirata 
(2004)M 
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